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ABSTRACT 

There is a maximum for the gravity of a black hole in the vertical direction in the accretion 
disc. Outflows may probably be driven from the disc if the radiation flux of the disc is greater 
than a critical value corresponding to the maximal vertical gravity. We find that outflows are 
driven by the radiation force from the disc if the dimensionless mass accretion rate at the outer 
radius m ou t > 1 (to = M /MrcM, M is the mass accretion rate, M^dd = ^Edd/O-lc 2 , and 
fvEdd is the Eddington luminosity). Assuming the outflow to be strong to carry away sufficient 
gas from the disc surface, we find that the radiation of the disc is limited by such outflows. 
The disc luminosity, Tdisc/^Edd oc In rh out , at large-m out cases. The Eddington ratio of the 
disc is ~ 3 for 77i ou t ~ 100, which is significantly lower than thatof a conventional slim disc 
without outflows (but it is comparable with that given in lKawaguchill2003h . This implies that 
the emission from some ultra-luminous X-ray sources with highly super Eddington luminosity 
should be Doppler beamed, or intermediate mass black holes are in these sources instead 
of stellar mass black holes. The spectra of the discs surrounding massive black holes with 
outflows are saturated in the high frequency end provided TO out > 2. We suggest that the 
saturated emission can be observed to estimate the masses of the black holes accreting at high 
rates, such as the narrow-line Seyfert galaxies, with the model calculations. The rate of the 
mass accreted by the black hole always M- m ~ M-Edd even if the mass accretion rate at the 
outer radius M out M-^dd, because most of the gas is removed into the outflows by the 
radiation force. If this is the case, the luminous quasars at high redshifts z > 6 should have 
grown up through persistent accretion at a rate close to the Eddington rate. 

Key words: accretion, accretion discs—black hole physics—quasars; general. 


1 INTRODUCTION 

Black hole accretion is the main energy source in different 
types of the objects, such as, X-ray binaries, and active galactic 
nuclei (AGN). In the past several decades, the measurements 
of masses of black holes in AGN with the reverberation map¬ 
ping met hod have been w idely used in the com munity (e.g., 


Petersonl 19931: Wandel. P eterson. & Malkanl 1 19991 : iKaspi et all 
200d : Peterson et al.l l2004h . The Eddington ratio distributions 


of large AGN samples show that the bolometric lumin osity of 
luminous AGN is roughly Eddington limited (e.g. jKollmeier et all 
120061) . or only a small fraction of AGN are super-Eddington 


(e.g.. Woo & Urrv 20021; Warner. Hamann 

& Dietrich 2004|_Wu 

2009 

; Shen& Kelly 2012|; Nobutaetal. 

20121: Kellv&Shen 

2013 

: Dong. Wu. & Cao 2014). but also 

see the discussion in 


Netzer & Trakhtenbrotl J20I4I) . The main observational features of 


AGN can be explained by accretion disc models quite well. The 
standard thin accretion disc model has been successfully used for 


the objects with moderate Eddington ratios, L/L^dd < 0.3 (L is 
the bolometric luminosity and L^dd is the Edd ington luminosity, 
IShakura & Sunvaevlll973l : lLaor & Netzeill 1 9891) . For more lumi¬ 
nous objects, the discs are no longer thin due to strong radiation 
pressure, and slim acc re tion d iscs are suggested to be in these ob¬ 
jects dAbramowicz et al.ll988l) . The slim disc models have been ex- 
tensiv ely stud ie d either analytically or w ith nu merical si mulati ons 


1999 

; Wang & Netzei 

20031: Gu&Lu 12007 

: Li. Yuan. & Cao 

2010 

; Sadowski et al 

201 ll: Eggum. Coroniti. & Katz 19881: 

Ohsugaetal. 2005; Ohsuga & Mineshige 200A 2011; Yang et al. 


the cooling timescale in slim discs dKawaguchU 120031) . A 
fraction of photons are unable to escape from the disc , and 
therefore they are trapped in the disc dBegelmaiii Il978l l. The 
radiation efficiencies of slim discs are lower than those for 
standard t hin discs d u e to the radial en ergy a d vection in the disc 

dSzuszkie wicz. M alkan. & Abramow iczll 19961 ; W atarai et alJhood : 
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IWatarail l2006lt . The detailed studies show that the luminosity 
of a slim disc increases slowly with mass accretion rate as, 
L/Lvdd oc lnm (m = M/Me dd, M is the mass accretion 
rate, and M^dd = T/Edd/O.lc 2 ), at high accretion rates, and the 
luminosity of the slim disc is alm ost saturated at around 10Z/Edd 
dWatarai et al]200Cil ; IWatarail2006h . 

Unlike most analytical works on slim discs, outflows are al¬ 
ways present i n the numerical simulations of super-criti c al ac- 
cretion flows JOhsuga et alj 1 20051 : lOhsuga & Mineshigel 1 20 lit 
I Yang et alJlioIiT This implies that only a sm all fraction of gas 
in the disc is finally swallowed by the black hole. lMcClintock et al.1 
d2006h found that the gas will be blown away by the radiation force 
if the relative half-thickness of the disc is H/R >1/v/2 (see equa¬ 
tion 7 and the discussion in their paper), in which an exact expres¬ 
sion of the vertical component of gravity is adopted instead of the 
approximation GMz/R 3 adopted in most literature, which is only 
valid in the case of H/R < 1. A series of works on the vertical 
structure of the slim disc indeed show that outflows are inevitably 
driven from the the disc surface provided the luminosity of the 
disc reaches the Eddingto n value dGu & "□ [20071 : iGu et alJEooi 


many luminous c 

uasars (e.g., Pounds et al. 2003]; Wang et al. 

2005 

Zhou et alj 2006 

; Gangulv & Brotherton 2008; Wang et al. 

2008 


20091: iTombesi et al.l20ld.l201 lL 120121) . 


In this work, we investigate the structure of super-critical ac¬ 
cretion disc with outflows driven by the radiation force of the disc. 
We describe our model calculations and results in Sects. [2] an dm 
The last section contains the discussion of the results. 


2 THE CONDITION FOR OUTFLOWS DRIVEN BY 
RADIATION FORCE 


It was suggested that the outflows are driven by the radiation 
force at the disc surface when the relative disc thickness H/R 
is gre ater than a certain value l lMcClintock et al.|l2006l : IGu & LtJ 
2007|). In this section, we extend the analysis on this issue in 
McClintoc k et al.l d2006ll . 

The equilibrium in the vertical direction of an accretion disc 

gives 

dp _ GMpz 

dz (R 2 + z 2 ) 3 / 2 ’ ’ 

at radius R. For a radiation pressure dominated accretion disc. 


where e r is the energy density of radiation. Substituting equation 
CD into 0. we have 

1 de r GMpz _ q(z)K T p 

3d 2 ~ “ {R 2 + z 2 ) 3 / 2 c ’ ( ’ 

where kt is the electron scattering opacity, and q is the radiation 
energy flux in 3-direction. We obtain 


, , GMz c 

q G) - {R2 +z 2 )3 /2— ’ 


(4) 


from equation 0. 

The half-thickness H of the disc is related to the outgoing 
radiation flux / ra d = q(H) from the unit surface area of the disc at 
z = H by 


/rad = q{H) 


GMH c 
(R 2 + H 2 ) 3 / 2 k t ’ 


(5) 


which can be re-written in dimensionless form. 


f , = - E -. 

(1 + H 2 ) 3/ 2 

The dimensionless quantities are defined as 


/rad — /rad 


R 2 K T 
GMc ’ 


and 


*-!• 


(V) 


The estimate of H in this work is different from and twice of H' = 
c s /Gk in most previous works of accretion discs (e.g., Frank, King 
& Raine 2002; Kato, Fukue & Mineshige 2008; see Section 4.1 for 
the details). 

The vertical component of gravity of the black hole increases 
from the mid-plane of the disc along z, and reaches the maximum 
at z/R = v/2/2 (see equation Q] and Fig.0. The balance between 
the radiation force and the vertical gravity is broken down if the 
radiation flux is sufficiently large. This implies that outflows may 
probably be driven from the disc if the radiation flux of the disc 
is greater than a critical value corresponding to the maximal 
vertical gravity, which is given by 


n max 
J rad 




(1 + HI ax ) 3 / 2 


2y/3 

9 


( 8 ) 


where ZZ max = \/2/2 is used. 

Both the vertical gravity and radiation force exerted on the gas 
in the disc is proportional to gas density p, and therefore the density 
p is canceled in the vertical force balance equation 0. We do not 
need to assume a specific function of z for gas density in the above 
analysis. This means our conclusion does n ot dep end on the ver¬ 
tical density distribution of the disc. iTurneil <120041) found that the 
energy dissipation is roughly constant vertically in the radiation- 
MHD simulation for a radiation-dominated accretion disc. If this is 
the case, the radiation energy flux in the disc is 

q(z) = / rad -2, (9) 

The radiation force exerted on the unit of mass in 3-direction is 

F(z) =q(z)EL = f reA l-9M-. ( 10 ) 

We plot the radiation force in Fig. |T| to compare with the verti¬ 
cal component of gravity. It is found that the radiation force is al¬ 
most in equilibrium with gravity in 3-direction for the low-/ ra d 
(i.e., low-ZZ) case. This_is the case considered in the radiation MHD 
simulation given in lTumei] < |2004|) . The vertical gravity —pfl^z is 
adopted in (Turner] d2Q04l Vs work, which is a good approximation 
for z/R <K 1. When 3/7? is sufficiently large (> 0.2), the vertical 
gravity deviates obviously from the linear relation of z (see equa¬ 
tion]]}. If the energy dissipation can still remain constant vertically 
in this case, the radiation force may be dominant over the vertical 
gravity at the upper lay of the disc (see Fig.0. It should be cau¬ 
tious on this conclusion, because the disc thickness adopted here is 
derived on the assumption of the vertical gravity —pQ/^z. This is¬ 
sue is to be resolved by the future numerical simulation adopting an 
exact function of vertical gravity. In the case of / ra d > /™the 
radiation force always dominates over the vertical gravity in the re¬ 
gion of z/R > Z7 max for arbitrary 3-dependent energy dissipation 
rate even if the disc thickness is not well determined. This can be 
verified by shifting / ra d (/rad > / r ™| x ) hr the horizontal direction 
in Fig.[T]for arbitrary value of H. 

The condition of outflows driven by radiation force is derived 
above for a radiation pressure dominated accretion disc, which is 
probably the case for a disc accreting at a high rate. We should 
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Limits on luminosity and accretion rate of a disc 3 




Figure 1. The vertical component of gravity as a function of z. The dashed 
line corresponds z/R = \/2/2. The coloured lines represent the radiation 
force in the case of the vertically constant energy dissipation rate with dif¬ 
ferent values of radiation flux, / ra d =0.5 (red), 2\/3/9 (green), and 0.1 
(blue). 


point out that the condition is valid even if the disc is not radia¬ 
tion pressure dominant. The condition for radiation pressure driv¬ 
ing outflows is the radiation force being greater than the maximum 
of the vertical gravity, i.e., 


,/rad P^T 


. -7-ir 

> pF z 


( 11 ) 


where T?_ is the vertical gravity for unit mass of gas, 
GMz 

•Fz ( 7 ?, z) — 


(R 2 +z 2 ) 3 / 2 ' 


( 12 ) 


of which the maximum is J r ™ ax = 2y/3GM/9R 2 . In dimension¬ 
less form, it becomes 


/rad > 2^3/9, (13) 

where equation 0 is used. This is exactly the same as the condition 
derived in the first part of this section (see equation 0. The self¬ 
gravity of the disc may dominate over the gravity of the black hole 
in the outer region of the disc, which is not considered in this work. 
This may alter the condition for outflows driven by radiation force, 
and the self-gravity may hamper the disc to launch outflows. 


3 THE ACCRETION DISC WITH OUTFLOWS 

For a disc accreting at a high rate, the outgoing radiation flux 
/rad > / r ™ d x - I n this case, a fraction of gas in the disc will 
be channeled into the outflows by the radiation force. The mass 
accretion rate in the disc will be a f unction of radius. Follow¬ 
ing the standard accretion d isc model IShakura & Sunvaevll 19731 : 
Frank, King, & Rairiel l2002l l. we derive the dissipation power in 
unit surface area of the disc with a 77-dependent mass accretion 
rate M(R), 


Figure 2. The radiation flux rates / rad from the unit area of the disc surface 
as functions of radius with different values of m = m C rit = 104 (red), 
0.2 (blue) and 2 (black). The red dashed line represents / rad = 


where A7i n is the mass accretion rate at the inner radius of the disc 
Ri n . For a non-rotating black hole, Ri n = 6GM/c 2 . The dimen¬ 
sionless form of equation 03 is 


D(r) = 

r 

where 

D{r) = D(R) 

7/Edd — 


1 _ min / r in\ 
rh(r) \ r J 


1 / 2 ' 


(15) 


R kt 
GMc' 

47 tGMc 


M 

Mvdd 


Mftdd = 


7/Edd 

0.1c 2 


Kt 


R_ 

R s 


R g — 


GM 


(16) 


In this work, we avoid being involved in the complexity of 
outflow physics. In the region of the radiation flux being greater 
than the maximal value given in equation 0, the radiation force 
dominates over the vertical gravity for the gas in the upper layer 
of the disc (z/R > T7 ma x), and it is blown away by the radiation 
force. This makes the mass accretion rate decrease, and therefore 
the energy dissipation rate decreases due to the decrease of mass 
accretion rate in the disc. We assume that this process continues 
till the radiation flux / rad — /™ d x - We neglect the radial energy 
advection, and assume /rad(-R) = D(R ) in the disc, which is to be 
justified later in this section. 

There is a critical mass accretion rate, m C rit, below which no 
outflow is driven by the radiation force. The critical rate is calcu¬ 
lated with 


D(r) 


15mcrit(r) 

r 



/•max 

Trad 


(17) 


at r = 9ri n /4, because D(r) reaches the maximal value at this 
radius for a given m. Using equation 0. we obtain 

merit = ~ 1.04. (18) 

5 

The radiation flux f ra ,d(R) varying with mass accretion rate rn is 
plotted in Fig. [2] 


D(R) 


3 GMM(R) 

\ Min (Rin\ 1/2 ' 

8tt R 3 

M(R) V R J 


(14) 


For an accretion disc accreting at a rate rh > m cr it, the out¬ 
flows are present in the region of the disc between r™ ln and r™ ax . 
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Figure 3. The outer radii of the region with outflows as functions of mass 
accretion rate m 0 ut- The solid line represents the result calculated with 
equation 02}- while the dashed line corresponds to that calculated with 
equation <H3. 


The gas at the disc surface starts to be removed into the outflows 
from the disc at r“ ax , where £>(r“ ax ) = /“jf 1 . The outflows are 
switched off at r™ m = 9n n /4, because D(r) will decrease within 
this radius (see Fig. 0. This means that the mass accretion rate at 
the inner radius of the disc mi n = m cr it if the mass accretion rate 
at the outer radius of the disc m ou t > m C rit- The outer radius r™ ax 
of the disc with outflows can be determined by solving the equation 
D(r 


15 m 0 


1 - 


TT^-out 


1 / 2 ' 


2y/3 
9 ’ 


with a given value of m ou t. Using equation GU- we have 

max ^ 45V^3 . 

^ w ~ T ^-out, 


(19) 


( 20 ) 


for the case of m ou t m cr ;t. The outer radii of the region with 

outflows varying with mass accretion rate m ou t calculated with 
equations GU and ( l20t are plotted in Fig. [3] It is found that the 
results almost converge when mass accretion rate m ou t > 5. 

The luminosity of an accretion disc with outflows (m > 
merit) can be calculated with 

oo oo „ 

I AnRD(R)dR = J ll^ldr, (21) 


T/Edd T/Edd 


where 


/rad (f) — < 


1 - 


(-) V2 ] ; 


r ^ C m : 

™ ax ; ( 22 ) 

r ^ r“ ax . 


r max 
J rad ? 

15m ou t N _ r h in / r t ,\1/ 2 1 
r [ rh Q ut V r / J ’ 

The integral in equation d2 1 b can be carried out analytically as 

1/2 


A = 


T/Edd 


35mi, 

27r in 


15m 0 ut 


1 - 


2mi, 


3m ou t V r. 


2%/3 / r; 


+ f 111 


(23) 


Figure 4. The relation of Eddington ratio A of the disc with the mass ac¬ 
cretion rate. The red solid line represent the result calculated with equation 
GU, while the red dashed line is the result calculated with the analytical 
approximation {24j. The black line indicates the results with standard disc 
model, i.e., the Eddington ratio varies with the mass accretion rate linearly. 
The green lin e represents the resu lt of the slim disc without outflows (see 
equation 19 in lWatarai et alj2000l) . The definition of m in their work is dif¬ 
ferent from ours, and their rh has been divided by ten to be consistent with 
our definition of rh. 


Substituting equation f20t into d23t , we obtain an approximation 
for A, 

2min 
3777-out 

, 2^3 5^3 . 

H-In —m ou t 

57%/3 8vT5 ( V 3 V /2 , , - 

= ~Tro - - - + —FT ln “5~ m ° ut ’ (24) 

162 675 \m ou t / 9 3 

in high-m out cases. We plot the Eddington ratio of the disc chang¬ 
ing with mass accretion rate m ou t in Fig. [4] in which the analytical 
approximation ( 124b agrees with the numerical result quite well if 
mout 1-5. It is found that the Eddington ratio A oc lnm out 
when m out is high. 

As discussed above, the outflows are present when accretion 
rate is higher than a critical value, which means that the mass ac¬ 
cretion rate is a function of radius. The outflows are driven from the 
disc region between r™ ln and r™ ax , and the mass accretion rate in 
this region can be calculated by 


35 . 2^3 

-min H- 

162 9 


/jK/3_ 
\ 45?77.ou1 


1 / 2 ' 


2\/3 . 5\/3 


D(r) = 


15m(r) 


1 - 


1 / 2 ' 


which reduces to 


m(r) = 


15 


-r + m 


min / l"in \ ' _ rmax 

• / T l j J rad ? 

m(r) \ r / 

. /r in \ 1 /2 2V3 , 9V2 

nin W J =l35 r+ — 


2%/3_ , 9\/2__ 1/2 


(25) 


(26) 


The mass accretion rate remains constant radially, i.e., m = m ou t 
in the region of r ^ r™ ax , and m = ih\ n = 3%/3/5 for 
r ^ r™ 111 = 9n n /4 = 27/2. In Fig.0 we plot the r-dependent 
mass accretion rate m(r) with different values of m ou t. The mass 
accretion rate in the inner edge of the discs remains constant, 
rhin = 1.04, due to the presence of outflows. 
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Figure 5. The mass accretion rates m(r) as functions of radius with differ¬ 
ent values of m 0 ut = 2, 10 and 100, respectively. The red line is the outer 
radius of the advection dominated region in the disc (see eg nation [27V 

IWatarai et alJ ( i200dV s calculations on the slim discs show that 
radial energy advection is important in the inner region of the discs, 
and the region of the disc with 

r < f a dv ^ 3m, (27) 

is advection dominated (see equation 16 in lWatarai et all200d . and 
note that their rh should be divided by ten to be consistent with our 
definition of rh). The advection dominated region disappears in the 
disc when rh < 2, which means that the advec tion is neglibile in 
this case (see equation 19 in lWatarai et al.ll200(Slt . In Fig. [5] we find 
that the mass accretion rate rh derived in this work is always much 
lower than that required for advection dominated case at any radius 
in the disc. This justifies the assumption of D(R) = / ra d adopted 
in our calculations. 

The effective temperature of the disc as a function of radius 
can be calculated with T S {R) = [f IS ,d(R) /cr] ly/4 by using the de¬ 
rived radiation flux (see equation [22}, when the black hole mass M 
is specified. The spectrum of the disc is therefore available with the 
derived effective temperature of the disc. We plot the effective tem¬ 
peratures of the discs surrounding massive black holes with differ¬ 
ent values of mass accretion rate m ou t in Fig. [6] The temperatures 
of the inner regions of the discs are almost the same for any values 
of m 0 ut > merit- In Fig. [7] we plot the radiation spectra of the 
accretion discs varying with mass accretion rates. The spectra of 
the discs surrounding massive black holes are saturated in UV/soft 
X-ray wavebands (v > 10 16 Hz) when accretion rate is as high as 
Wont ~ 2. 

4 SUMMARY AND DISCUSSION 
4.1 Structure of the disc 

The vertical component gravity increases with 2, and reaches the 
maximum at z/R = v/2/2 (see Fig. [TJ. The radiation force will 
dominate over the vertical gravity in the disc if the radiation flux 
is sufficiently large. It implies that outflows can be driven from 
the disc if the radiation flux of the disc is greater than a criti- 



Figure 6. The effective temperatures of the discs as functions of radius 
with different values of m 0 ut = 0.2 (black), 1 (red), 2 (green). 10 (yellow) 
and 100 (blue), respectively. The solid lines indicates the results calculated 
with black hole mass M = 10®Mg, while the dashed lines are for M = 
10 9 Mg. 



Figure 7. The spectra of the discs surrounding a black hole as functions of 
radius with different values of m 0 ut = 0.2 (black), 1 (red), 2 (green), 10 
(yellow) and 100 (blue), respectively. The solid lines indicates the results 
calculated with black hole mass M = 10® Mg, while the dashed lines are 
for M = 10 9 Mg. 

cal value. This has already been discus sed in the previous works 
dMcClintock et al]2006l:lGu & Lull2007l) . 

The relative disc thickness of H/R corresponds to a dimen¬ 
sionless radiative flux / ra a at the disc surface for a radiation pres¬ 
sure dominated disc. We derive the critical mass accretion rate 
riicrit = 3\/3/5 — 1.04, below which the relative thickness H/R 
of the disc is always less than \/2/2, i.e., no outflow is present 
(see Fig. [^. In the case of m > m cr it, the outflows are driven 
from the disc surface in the disc region between r™ ln and r™ ax , 
where r™ ln = 9n n /4, and r™ ax is a function of mass accretion 
rate rh (see equation [T9}. The outer radius of the disc region with 
outflows r“ ax oc rh when m is large (see equation [20}. For a disc 
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accreting at an extreme rate m ~ 100, the disc region with out¬ 
flows can extend to thousands of gravitational radii (see Fig. [3}. 
Due to the presence of outflows, the mass accretion rate is not a 
constant radially in the disc. A fraction of gas at the disc surface is 
blown away by the radiation force, so we assume that mass accre¬ 
tion rate is self-adjusted to let the dissipation power in unit surface 
area of the disc D(R) = /™ d x . This means that we have assumed 
the radial inflow timescale at a certain radius to be longer than the 
timescale for an outflow to change the accretion rate of the whole 
disc. The outflow timescale can be estimated only if the dynam¬ 
ics of the outflow is well studied, which is beyond the scope of 
this work. The mass accretion rate as a function of radius is de¬ 
rived in this work (see equation[26] and Fig. |5j. At large radii, the 
m ass accret i on rat e rn oc r, which is the same as the result given 
in lGu & Lul d2007l l. This is qualitatively consistent with the hydro- 
dynamic al sim ula tions on su per-critical accretion flows (see, e.g., 
Fig. 6 in lOhsuga et al.ll2003) . Their numerical simulations indicate 
that the outflow can be so strong to carried sufficient gas from the 
disc, though the detailed study of the outflow dynamics is needed 
to clarify this issue. 

We note that the relative disc thickness H/R is larger than 
that deri ved in the p reviou s works for the same mass accretion rate 
m (e.g.. IKawaguchil [20031) . In most of the previous works, an ap¬ 
proximation of the vertical component of gravity, zfl^, is adopted, 
which over-estimates the vertical gravity compared with the exact 
one u sed in this work when z/R is large (see Fig. 1 in iGu & Ltj 
12007b . This approximation leads to underestimate of H/R in the 
previous works. For relative disc thickness H/R < 0.2, the ver¬ 
tical gravity zQ.^ is a good approximation (see Fig. [TJ. Our cal¬ 
culations show that the disc th i cknes s H is systematically larger 
than that derived in lKawaguchil I2OO3I) . The relative disc thickness 
H/R —» V2/2 in the disc with R < 400 R e when m = 10 (see 
Fig.0, which is about twice of that given in lKawaguchil d2003h for 
the same accretion rate (see Fig. 5 in that paper). In that work, the 
scale height of the disc is estimated with 

H'=£-, (28) 


where the sound speed at the mid-plane of the disc c s = 
(pmid/pmid) 1 ^ 2 , which is different from ours even in the limit 
of H -C R- In this work, the disc thickness can be derived from 
equation 0. 


H ~ 


/radKT 

^K C 


when H/R -C 1. The outgoing radiation flux is 

_ 4qcT* id 

J rad n r r 

OKTPmidT/ 

Substitute equation ( 130b into ( 129b . we have 


(29) 


(30) 


H ~ 


(4p mid ) 1/2 

Pmid^K 



(3D 


where p m i d = aT ^ id /3 is used for a radiation pressure dominated 
disc. This can exp lain the systematic difference of disc thickness 
between fKawaguchi] d2003l) ’s paper and this work. 

Our calculations show that outflows can be driven from 
the outer region of the disc, which may be truncated at _B ou t 
where the disc b e comes gravitatio nally unstable (e.g., iHurell 19981: 
IKawaguchi, Pierens, & Hura 120041 ). The outer radius of the disc, 
Rout, is a function of black hole mass, accretion rate, and the vis¬ 
cosity parameter a, which can be calculated with derived disc struc¬ 
ture. However, the detailed disc structure has not been derived in 


this work, and the constraint on the size of a disc with outflows 
will be reported in our future work. Our present work cannot tell 
whether the outflow is damped by the vertical self-gravity or not. 


4.2 Luminosity of the disc 

The luminosity of the disc with outflows can be derived with equa¬ 
tion ( 123b . We find that, L/L^dd oc lnm, if m is l arge (equation 
24}. This looks simi lar to the result for the slim disc l lWatarai et al.1 
200(3; [Watarai 200(3), however, the physics is completely different. 
In their works, no outflow is considered, and the mass accretion 
rate remains constant radially in the disc. 

The cooling timescale is longer than the viscous timescale of 
the slim disc. Therefore, only a fraction of the dissipated power is 
radiated away, and the remainder is advected inwards. The radia¬ 
tion of the slim disc is limited by the photon trapping effect within 
a certain radius, which is roughly proportional to mass accretion 
rate m. Thus, the luminosity of t he slim disc L/L^dd oc lnm 
d Watarai et a l. 2000; Watarai 200(3). In our present work, the mass 
accretion rate at the inner radius of the disc rh; n = m cr ;t ~ 1.04 if 
m ou t > m cr it, because most of the gas in the disc is carried away 
in the outflows. The luminosity of the disc is therefore limited by 
the outflows (see Sect. 0. We find that the luminosity of the disc 
with outflows L /Z/Edd increases more slowly than that for the slim 
disc (see Fig. 0. In this work, the disc luminosity L/L^dd ~ 3 
provided m ~ 100, whic h is lowe r than L/L^dd ~ 10 for the 
same accretion rate in I Watarai et al.l d2000b . Our results are compa¬ 
rable w ith L/ Z/Edd ~ 2.6 for m ~ 100, as given in iKawaguchil 
d2003h . 

Our calculations are carried out based on the balance between 
the radiation force and the gravity in the vertical direction of the 
disc, which means that the main conclusions on the maximal disc 
luminosity in this work will not be altered whether the energy ad- 
vection is important or not. We find that the radial energy advection 
is always negligible in an accretion disc with outflows (see Fig. [5] 
and the detailed discussion in Sect.0. 

The UV/soft X-ray spectrum of the disc surrounding a mas¬ 
sive black hole is determined by the effec tive temperature profile 
and electron scattering effects dKawaguchill2003l) . In this work, the 
spectral calculations are carried out based on the local blackbody 
assumption without including the electron scattering effects. There¬ 
fore, the spectral shape of the disc in the high frequency end is pre¬ 
dominantly determined by the structure of the inner region of the 
disc. In Fig. [5] we find that the structure of inner region (< 100Z7 g ) 
of the disc remains unchanged when m ou t > 2, from which most 
of the accretion power is radiated. The effective temperatures of 
the discs with outflows in the inner regions almost converge for 
different values of m ou t (see Fig. [6}. Therefore, it is found that the 
emission in the high frequency end of the spectrum (v > 10 16 Hz) 
is saturated when rh ou t > 2 (see Fig.0. The further calculations 
of the spectrum including the electron scattering effects, as done by 
iKawaguchil d2003h . are necessary for detailed astrophysical appli¬ 
cations. As the structure of the inner region of the disc with out¬ 
flows accreting at a high rate is insensitive to the mass accretion 
rate at the outer radius, we believe that the feature of the emission 
saturation should still be present, though the final results may be 
quantitatively different from our present results. Such calculations 
are beyond the scope of this work, which will be reported in our 
future work. 
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4.3 Implications of the model 


The saturation of the emission in the high frequency end of the disc 
spectra can be used for estimating the masses of the black holes 
accreting at high rates, such as, narrow-line Seyfert galaxies. It im¬ 
plies that luminous quasars with m out > 2 can be used as a type of 
“cosmological candles”, and one does not need t o search mass ive 
black holes accreting at extremely high rates JWang et al .1120 1 3l) . 

The relation of m ou t-T//TEdd derived in this work shows 
that the radiation of the disc can only be several times of Edding¬ 
ton luminosity even if the mass accretion rate m ou t )§> 1. This 
seems to be in contradiction with the observations of some ultra- 
luminous X-ray (ULX) sources, of which the luminosity can be 
highly super Eddi ngton if normal stellar mass black holes res ide in 
the U LX (e.g., Feng & Sorialboi ll : ISutton. Roberts, & Middletonl 
120131: ISo ria et a l]|2014j , and the references therein). As suggested 
by King et alT i 200 ll) . this issue can be alleviated if the X-ray 
emission from the ULXs is mildly beamed to the observer, and 
the obs erved X -ray emission is predominantly from the outflows 
feegelman. King. & Prinel3l2006l) . which are qualitatively consis¬ 
tent with our results that outflows are inevitably driven from the 
disc surface when rh ou t > merit- The detailed investigation on 
such outflows from the radiation dominated accretion disc will be 
reported in our future work, which is beyond the scope of this pa¬ 
per. An alternative possibility is that these ULX may have interme¬ 
diate mass black holes, of which the lu minosity is therefore lower 
or close to Eddington luminosity (see lFeng & Soriall201 ll . for a re¬ 
view, and references therein). 

Luminous quasars are discovered at redshifts z > 6 (see 
Fan et all 120031 : iMortlock et al.l 1201 ll : IVenemans et alj 1201 3l . and 
the references therein), corresponding to the age of the universe 
less than 1 Gyr. To grow a black hole to ~ 10 9 Mq in such a 
short period of time, the massive black holes in these quasars must 
have either accreted persistently at a rate close to the Eddington 
value te.g.- lShapirdbopj ; [Tanaka & Haimanll2009t> . or unde rgone 
many episode super-Eddington accretion events (see King[ 2003; 
Kawaguchi et al.l 20041 : IVolonteri & Reed 120051 : IVolonteri & Silkl 


2014, nd the references therein). In this work, we find that the rate 
of mass accreted by the black hole min ~ 1, even if the mass accre¬ 
tion rate at the outer radius of the disc m ou t 1 , because most of 
the gas in the disc is removed in the outflows from the disc. Our re¬ 
sults seems to be consistent with the continuous acc retion sce nario 
for massive b lack hole g rowth at high redshifts (e.g., I Shapiro! 20051 : 
iTanaka & Haimanll2009l : lTanakall2014l) . 
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